


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1976 


An optimization study of a low thermal 
potential power system 


Buckingham, James Robert 


Monterey, California. Naval Postgraduate School 


http://ndl.handle.net/10945/17886 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
D U DLEY research materials and institutional publications created by the NPS community. 
get Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS'‘s first 
KNOX appointed — and published — scholarly author. 


LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 





_ 
= 
oc 
i Ss, » is 
2) RS 
ad 
b Saez ean 
Les 
<= 
= 
oe 
a. 
& 





-QUDLEY KNOX Lim RARY 





NAVAL POSTGRADUATE SCHOOL 


Monterey, California 





THESIS 


AN OPTIMIZATION STUDY OF A LOW 
THERMAL POTENTIAL POWER SYSTEM 


by 


James Robert Buckingham 


june 1976 


Thesis Advisors: M. D. Kelleher 
Wm. M. Raike 





Approved for public release; distribution unlimited. 


T1/4018 








SECURITY CLASSIFICATION OF THIS PAGE (When Date Entered) 


REPORT DOCUMENTATION PAGE 


- REPORT NUMBER 2. GOVT ACCESSION NO 


4. TITLE (and Subtttie) 


n Optimiz zation Study of a Low Thermal| Master's Thesis 
otential Power System June 1976 
6. PERFORMING ORG. REPORT NUMBER 





READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


5. TYPE OF REPORT & PERIOD COVERED 























7. AUTHOR(e) 8. CONTRACT OR GRANT NUMBER(e) 







Ledr. James Robert Buckingham 





10. PROGRAM ELEMENT, PROJECT, TASK 
AREA & WORK UNIT NUMBERS 


12. REPORT a 
June 19 
13. tO. OF — 


18. SECURITY CLASS. (ol thie r@port) 
Unclassified 





$. PERFORMING ORGANIZATION NAME ANO ADORESS 


Naval Postgraduate School 
Monterey, California 93940 













- CONTROLLING OF FICE NAME AND AOORESS 


Naval Postgraduate School 
Monterey, California 93940 





from Controlling Office) 


Naval Postgraduate School 
Monterey, California 93940 
















Se. DECL ‘oulte OOWNGRADING 


SCHE 
- DISTRIBUTION STATEMENT (of thie Report) 


** Approved for public release; distribution unlimited. 


- DISTRIBUTION STATEMENT (of the abstract entered in Block 20, if different from Report) 


- SUPPLEMENTARY NOTES 





»- KEY WORDS (Continue on reverse side if neceesary and identify by biock number) 
Nonlinear programming, thermodynamic cycle, ocean thermal 
energy conversion, system optimization, cost estimation, 

system engineering. 












. ABSTRACT (Continue on reveree side if neceseary and identify by bleck number) 
A power generating system using the low thermal potential 
available from the vertical temperature distribution of the 
ocean is analyzed as a combined engineering and economic 
mathematical model. The model is optimized for minimum 
capital cost employing a sequential unconstrained minimization 






DD , ANTS 1473s EOI TION OF 1 NOV 68 1S OBSOLETE 
(Page 1) S/N 0102-014 6601 | an 
SECURITY CLASSIFICATION OF THIS PAGE (When Date Entered) 





a NSN (ASR gery pA YE SESE 
Se CURITY CLASSIFICATION OF THIS PAGE(Whon Dota Entered. 





algorithm. Examples of the kinds of engineering and cost 
information available from the model are presented. 


DD Form_ 1473 


an : 
S/N 0102-014-6601 SECURITY CLASSIFICATION OF THIS PAGE(When Data Entered) 





AN CETIXZIZATION STUDY ee cel" THERMAL EFOTENTIAL EQWER 


by 


James Robert Buckingham 
Lieutenant-Commander,”’United States Navy 
B.S., University of Illinois, 1966 


Sukmitted in partial fulfillment cf the 
reguirements fer the degrees of 


MASTER OF SCIENCE IN oe ENGINEERING 
an 
MASTER OF SCIENCE IN MANAGEMENT 


from the 
NAVAL POSTGRADUATE SCHOCL 
June 1976 


“Thesis 
BIas 
c.f 





poe KNOX LIBRARY 
AVAL POSTGRADUATE 5 
CHOOL 
MONTEREY. CALIF. 93940 


AESTRACT 


A pewer generating system using the low thergal 
frotential available from the vertical temperature 
distrikuticn of the ocean is analyzed as a combined 
engineerirg and economic mathematical model. The 
model is cptimized for minimum capital ccst employing 
a seguential unconstrained minimizaticn algorithn. 
Examples CE the kinds of engineering and ccst 
inforgaticn available from the model are presented. 





TABLE OF CONTENTS 


emer leE ee aU eet Ue sien a. ste daisliniwiiajlccaileilells( 6. elisliellieievepeiee ses 600 0 ee dy 

cd ERMC OEMOMNLiclelcls swe ce «closes ccocscsscocssesceceueece 2 

isle el (GIS rdf 0e |, iat Sh aea eae aa ire inane ia a ar ee ae 

Dee ee EOLA) O ORCrEo LU uelnMalsietstolcls 4 s 6 » s\sisliéieie < ejleldle «0 os < sees 6s “ab 
Mo GRR DOOR LEON sc ce wee ec cect sees anes a ees as Zils 

Eo Gee DNA tsetse we cc ce cc cee sae e cece een eee ee Zo 

C. ENGINEEFING FRAMEWORK... cc cence ccc cece cece ee § 30 

1. The Borler..cccccccccccceveccccccesececces 30 

7a gis) (SOCAN SE OTC tite SCO eee I nS 

Se SS PbS tereeie a) 6 < = s seieleess 2s sees ees tse ene ee 

Tie CEST PNRPREMNINC CMM ese 55 in5 15:5,0).59aslvuana is leiwiaie se Ss 4y4 a) LEO 

ES ASSEBELY OFT THE valnenAleecaAlL MODEL. SC. . ses 4s es 53 

1. Objective Function and the Censtraints... 55 

ee the optimization MEthod.......4. se eees se Neen 

Oem. RESMRIS RIND COMERUIGINOINS . 0.5... cs sce ewes ese cee «OZ 
ee Ger eueeete he © Neo etc s oleielniel oeletol ela <)e!elcl alse «s/s 6 24 «)s =e" ° val 
Appendix 2: CCMPUTER PROGRAM AND SAMPLE DATA.......ce. FAG 
ieeeresen DERE PRNGES .. so. Scleugeiee sees ccc sec os 550s en ce co cle wae MOS 
Teen CISITEIEUTICN LISD. sone osc cess c cece cctseeseceees | Lee 





Eist Of TABEES 


Material Factcrs for the Heat Exchangers..... 
Material Factors for the Centrifugal Pups... 
Material Factors for the Propeller Pumps..... 
Defiriticn of the Elements of the Vectcr X... 
Prom SRMi np Deas cc ese ce ct ee te Rete cae es ee 
Example Exoblem Solution Summary. ....ccccecee 
SOmucacn POWET BMAGEt CS .. 2. os cs cc 0 oss oe ee sees 


Additicral Data’ from Sample Problem Sclution. 


pill 
51 
52 
54 
64 
65 
68 


10 





LIST OF FicuRes 


Some @ Gc vem ti C DEAGEAM << .6 cers oo 6s « 0) 0-0 os tenaianeraie 6 o's 
Sp Slecme f lela CWiciiiiietiel elo «0 so s claisle os «6 «2 0 6s clsieiels cos «ss 
CPA COV G URI GiVGHEG os 6s a ale 6 s\sisis oo 6s ss os 6 secs s sie slates - 
Biermeer Ch@Gi= SNOT 1ON.6. s.0 1s cleo <c.0 ces siecle eis oe) o eels «ais 


Compariscn of a Single Tube Versus 


a Tucke in a Bundle. @eeeeenes#?steseeas@se#ss#tee@eeeee#eee*e@egeeee@eeee#ee@eese#s8eee?ee 
Thergzal Network for a Single. @®eg@eee@e%e?ed@seeoeodeoe@eeede?e@e#e@ee#eee#e?8 e820 @ 
Pmomere@ SCIULION SEGUENCE.. ..cass cewe cceenaac seas sels 


Cost EL eGaediow rn sc. « Seetees SU. 0 0 0 cee Ge Ce wee se ce eels sce eee 


16 
Ze 
27 
Sie 


oO” 
De, 
Dy, 
67 





axe cOLS 


in 


© oo ee oe 


cm ae Bb we 


hs 


Kc 


Ke 


Nsg 
NTU 


LIST OF SYMEOLS 


CDESCRIPTICN 


area 

heat capacity 

C@apacity rate 

tuke diameter 

diameter of the shell 
Fanning friction factor 
factcr for the tube 
bundle boiling 

Inequality constraints 
equality constraints 

heat transfer coefficient 
erthalry 

erthalrpy corrected for 
sukcccling 

height frem midplane of 
tke Eciler to top of 

the tuke kank 

vapor quality 

thergal cenductivity 
pressure loss factor at 
the tuke 
pressure 
tke tuke exit 
tuke length 
mass flow rate 


entrance 
loss factor at 


averace numter of tubes in 
a ccivmn in the condenser 
Specific speed of a guar 
nugber cf heat transfer 
units cf an exchanger 
Nusselt numker 


TLECSsvce 


UNITS 
ENGLISH MEDRIC 
ine cné 
BTU/lbm-°F J/kg- °C 
BTU/hr-°F kW/°C 
in cn 
fete m 


BTU /f£te—-hr-PF 


BTUy/ilbo 
BTUslbo 


oe 


BLU/RE-ft e+ 


Et 
lbmyhr 


RPM 


lbtyin2 


kW/me-9C 
J/kg 
J/kg 


kW/a-°C 


kg/hr 


RP 





Fi 
EL 


ke 


A Wm tA 
Do 


a G&G wet 


VL 


GREEK 
SYBEOLS 
E 


? 
A 


© Qo < © 


PpCwer 

tue ritch 

Prandtl numtker 

heat flux 

heat transfer rate 

heat transfer 
(R1,R2,R3,R4) 


Reyncids nugber 


resistance 


entrory 

number cf shells 
segmental area 
ture 
terperature 


wall thickness 


velocity 

cverall heat transfer 
ccefficient 

specific vclume 

vapcor lcad 
propertionality factor in 
the transition region 
capital cost 


heat transfer errectiveness 
efficiency 

Prepertionality factor in 
the transition region 
dyraric viscosity 

kinematic visccsity 

density 

surface tensicn 

free flew/frontal area 


ratio 


BTUyshr 

rial 

BTU /hr-ft2 
BTU/hr 
hro-°F /BIU 


BTU/Lbm-98 
fte 

in 

OF 

ft/sec 


Devitt 2=npoor 


ftsyibdn 
ibm/hr-f£t43 


bed 


dollars 


lbm/sec-ft 
ft-ysec 
lbmyit3 
Bene ays 


MW 


cn 


kW/me 
Mw 
OC/KW 


kW-hroykg-°K 
m< 

cn 

°C 

m/sec 
kWysm2-9C 


m3/xg 
KG7 aE ale 


doliars 


N/Ssec-m 
me/sec 
Kg/m3 
N/® 





SWESCRIPTS 


ab 
E 
bulk 


cal 


cp 
cyc 


fg 


12a 


he 


sal 


cut 


req 
Si 
sys 
LE 
Tu 


aksorked 

bEciler 

ayerace of the inlet and 
Outlet temperatures 
ccndeérser 

critical pressure 

the cclder fluid in 

the heat exchanger 

carnct efficiency 
circulation puap 
therascdynanic cycle 
Saturated liquid 

change in the guantity as 
tke fluid changes from all 
liguic to all vapor 

feed pfunp 

Saturated vapor 

the hetter fluid in 

the heat exchanger 

heat exchanger 

inside 

fluid entering the excnanger 
laginar flcw regime 

mean temperature difference 
max. reossitle temperature 
dirference 

outsice surface 

fluid leaving the exchanger 
at ccorstant pressure 
rejected 

seawater 

the system 

transition regime 


turkulent regimse 


10 





tur turkine 
W wall 


wf acrking fluid 


¢ maximtn temperature difference 
1 State point 1 

zZ tet ae Dy 

zs . ” 2s 

3 at ae 

4 m " 4 

5 1 1 5 

ss ” : 5S 


11 





A. EACKCECURNT 


As cf 1874, the United States contained only 6% cf the 
world pogulaticn but was using 33% of the energy consumed 
€ach year according to Ref. 1. Approximately one third of 
the oil used in the U.S. is imported. This heavy dependence 
cn foreign scurces fer energy coupled with the exhausticn of 
dcrestic fessile fuels, particularly oil, in the foreseeable 
future is causing the U.S. military and civilian sectcrs of 
the economy tc search for energy sources that would be 
inexhaustible and independent of fcreign ccntrol. 
Increasing fuel costs and environmental prcklems are other 
forces fushing the search for alternatives. Nuclear fower, 
once considered the cure-all, is facing severe frotlems. 
Ccests cf ccnstruction, operation, and fuel have risen 
dramatically. Questicns concerning the safety of cperation 
of the nuclear plants and cf the storage cf nuclear sxastes 
have beccrwe felitical issues that have slowed the 
exploitaticn cf nuclear power. The present nuclear reactors 
use a. fuel tkat, like petroleun, has a limited availakility 
and whose price is increasing rapidly. The breeder reactor, 
which would create more fuel than it burns, and the fusion 
reactor, which would use the hvdrogen found in water, are 
still a lcng way frem producing power for a world whose 


demands fcr energy increase each year. 


Many fecfle advocate the explcitation cf the sc-called 


"free" enrercy scurces. These scurces are: (1) the kinetic 


V2 





energy of mcving fluids, such as winds and ccean currents; 
(2) the fetential energy of tides and rivers; (3) the heat 
generated inside the earth; (4) the direct conversicn of 
sclar energy into electricity and heat; and (5) the use of 
Sun warmed water in conjunction with a scurce of cclder 
water te prcvide the temperature differential to run a vapor 
power cycle. Nene cf these sources are truly "free". What 
is happeninc is the trading off cf the transportation, 
prcecessing, and environmental costs of a concentrated energy 
scurce fcr the capital, ofgerating, and social costs of 
convertinc a very diffuse source into a mcre concentrated 
fcrm such as electricity. 


The earth can ke thought of as a giant heat engine, 
abscrbing energy frem the sun and reflecting energy as 
thermal radiaticn kack to space as fpictured in Ref. 2. 
Since the feles receive less energy than the equator, the 
atmosphere and ocean attempt tc distribute the energy more 
evenly cver the earth. The air is heated Fy absorption of 
sunlight and Ey contact with the surface of the earth. The 
heated air rises thus setting up surface wind currents as 
ccecler air tries tec replace the rising air. The shear 
stress created by the relative moticn between the air and 
the surface water causes the water to move thus creating the 
Surface water currents. Many currents in the oceans, 
recardless cf depth, originate with the wind shear stress at 
the air-surface interface. The circulaticn patterns cf tne 
oceans beccme extremely complex due to the influence cf the 


rotation cf the earth and the shapes of the ccean basins. 


In general the surface waters near the equatcr are 
wargzed ty the sun and flcw towards the poles giving uf 
energy by é€vapceration, and by convecticn tc the air ard by 
radiaticn tack to space. AS these currents cool, scmne of 
the water teccmes dense enough tc sink and then flows back 


to the eguatcr alcng the bottcn. 
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The relaticnship between salinity and temperature and 
the density cf seawater is such as tc create a stable 
stratificaticn of the ocean in many areas. A salinity 
decrease cr a temperature increase reduces the density. The 
tetal effect is to create a warm surface current overlaying 
the cold water returning from the poles. In the trcepical 
waters the surface temperature may exceed 80°F (2€.7°C) 
while 3=,00C to 6,000 ft. (914 to 1828 m) dcewr the 
temperature may be about 40°9F {4.4°C). This temperature 
difference can te used to rur a man-made heat engine, like 
the Ccean Thermal Energy Converion (OTEC) power plant 
ccrcept. 

The proper choice of a site is critical to the ccst of 
such a heat engine. Since the temperature difference is low 
Ccupared to ccnventional thermal power generating methods 
which typically use temperature differentials exceeding 
50C°F (2€0°C), the effect of a 19°F (0.56°9C) loss in the AT 
is much nore costly at the low thermal differentials, which 
are usually less than 509F (27.8°C) in the cc¢can. 
Therefore, the vertical temperature profile of the cceans 
will ccntrel the location of an OTEC power plant. The 
Surface temperatures vary with location and with time. At 
mcst places in the ocean there is a seasonal change of the 
surface water temperature of 5 to 89°F (2.8 to 4.4°C) except 
in the trepical areas that remain relatively corstant 
acccrding to Ref. 3. However, ali regions cf the oceans are 
affected Ly wind-yenerated mixing of the surface water that 
Changes tke surface temperature by a degree or sc anda 
daily cycle cf a few tenths cf a degree. Some cf the 
pessibie U.S. near-shore sites being considered are the 
waters arcund the Bawaian Islands and the Gulf Streaz off 


the scutkeastern U.S. 


In the past, the engineering analysis and the cost 


estimaticn tave Ereen conducted as if they were separate 
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furctions. In a systen With a stall temperature 
differential, such as this one, the engineering desicno and 
eccncmic analysis should be iinked tcgether as a single 
precedure because of the capital intensive nature cf the 
ocean therral power plant. For example, the heat transfer 
egtaticn (¢€ = UA AT) is coupled to the ccst eguaticn for 
heat exchangers ($ = K A") through the surface area of the 
tukes. Eecause cf the small AT involved the area mtst be 
very larce tc supply the same Q. The other components are 
all linked tccether by the engineering and ccst equaticns so 
that the preklem cannot logicaily te fEroken intc an 
equipment cy equipment optimization. The present analysis 
considers a clesed vapor power system oferating cn the 
thermal fetential available from the ocean. The eccnomic 
Ofptimizaticn cf the system is carried out to shew what 
informaticn can te gotten, what conclusicns can be reached 
frcm suck analysis and where research effort shculd be 
expended tc imprcve the design of the system. This study 
does not attenzrt tc predict the cost of such a system cr to 


make specific recosumendations about the parareters. 


A schematic diagram of the system is shcwn in Pigure bie 
The warn surface water, the Gulf Stream for example, 1s 
pumped thrcuch a heat exchanger. There scme enercy is 
transferred Ey boiling the working fluid such as amzonia. 
The working fluid is piped to a turbine where some of the 
energy is ccnyerted to mechanical energy then to electrical 
energy in a generator. The unavailable ortion of the 
akscrbed energy is rejected to the environment in the 
condenser anc the working fluid is returned to the Mliguid 
state. The lcw temperature in the condenser is maintained by 
EFueEing ccld ocean water up from the depths. Finally the 
working fluid is pumped back to the boiler. 


The Ccear Thermal Energy Ccnversion program has received 


much attenticn in the last several years from many 
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researchers funded at first primarily by the Naticnal 
Science Fcurdation (NSF). Since the estaklishment cf the 
Energy Research and Development Administraticn (ERDA), new 
research ccrtracts are funded mainly by this agency. The 
C.S. Navy is alsc funding some research in the OTFC area. 
Ref. 4, 5S and 6 give a good summary of the work that has 
been done. thile much of the work has concentrated on the 
heat exchangers, the other components and the system as a 
whcle are receiving attention. Some researchers are 
concentratin¢c on the environmental and legal gquesticns of 
such systems. Cther teams are considering hcw to use the 
energy supplied from an OTEC fower plant. 


In any désign, trade-offs must be made. Capital cost 
can be traded for operational costs. Pumping power can be 
traded fcr heat exchanger size. Trade-coffs can involve 
items that kave a common measure such as dollars; scme other 
trade-offs are very hard to measure in commen terms ‘such as 
sccial and environmental costs. The trade-cffs may involve 
crly one ccmzpcnent such as a pump or a turbine, but most 
Often the trade-cffs affect more than ocne compenent or 
pessikly tke entire systen. 


The trade-cf£f process shculd continue until an cretinun 
design is crtained. From Ref. 7, "The aim cf optimization 
is the selection, cut of the multiplicity of fotential 
sclutions, of that solution which is the best with respect 
to some well-defined criterion." In crder for the plant 
dé€sign tc be cptimum, the compenents must be cptimized 
within the ccntext of the entire syste. Some cf the 
parameters fcr a piece of equipment may have no ¢cirect 
effect cn the cptimum design of another porticn cf the 
System, Eut ¢cme parameters will have significant effects. 
For exagple, a heat exchanger is made up cf a large number 
of tubes cf a given diameter and length and is designed to 


have a certain heat transfer rate for a given flow rate of 


I/ 





water pumped through the tubes. The diameters, length and 
the numker cf tubes may be optimized for lowest cost cf the 
heat exchanger. This approach leaves out the effect that 
the tube diameter and the tube length have cn the pressure 
drcp acrcss the heat exchanger. Minimizing the cost cf the 
teat exckarger irvolives minimizing the surface area. This 
may lead tc increasing the pressure drcep acrcss the 
exchanger. The higher the pressure drop the more forerful 
and the gecre expensive the pump must be to maintain the flow 
velccity. This is just one example cf the hundreds of 
trade-offs fessible. 


Unfortunately, what happens in many cases is that each 
grcup cf experts werking on a particular ccmfonent optimize 
that iter arcund some given conditions that may have Leen 
set ky teams working cn some other section cf the systen. 
Little analysis is made of how the entire system resrecrds to 
the coupling ccnstraints. All that the purp manufacturer 
wants is that the bLuyer tell him what the flow rate, head 
and service ccnditions are and he will quote ae frice. 
Befcre tke buyer can decide on the flow rate and head he 
sheculd first know how the ccst cf the pump varies as flow 
rate and tead vary so that those two parameters can be 
optimized in tke ccntext of the systen. 


The epretlem then becomes cne of acquiring sufficient 
infcrmaticn cn hcw costs vary as certain parameters change. 
Each manufacturer is able to make estimates cf costs fcr his 
prcduct but te is reluctant to devulge his informaticn to 
kuyeé€rs cr ccempetiters. This makes constructing cost curves 
a hazardcus precess at best, since the data nust be gathered 


by cther mwetkcds that may cause large unexplained variances. 
Sometimes the optimizaticn analysis 1s made with 


assumpticnos that make the answers of guestionakle value. 
For example, sometimes the heat transfer coefficient, U, is 


ce 





assumed fixed in the equation g = UA AT. But, as shcwn in 
Ref. &@, for toiling the heat transfer coefficient is very 
sensitive tc the temperature difference with the result that 
gq = UA'TAT3S- 33. The The accuracy of the engineering eguations 
shculd re kept in mind. Most of the equations are 
ccrrelaticns tc fit experimantal data. The correlation was 
develcped fcr scme particular set of data and may disaaree 
Significantly when ccmpared to data taken by cther 
experimenters. The limitations of the various relations 
must be kept in mind when reviewing the results. For a 
system the size of the OTEC plant, a pilct plant should be 
ccnstructed in order to prove out the answers given by the 
calculaticns refore expensive mistakes are made. 


As a cesign progresses, the analysis begins with a 
general ficttre ccntaining many Simplifying assumpticrs and 
prcceeds foreward with more details added at each stage. 
This report is the second one written at the Naval 
Pestraduate fSchccl on the subject of OTEC. Commander Furman 
Sheppard, USN, in his analysis, Ref. 9, attempted tc show 
what kinds cf information could te gained from a corbined 
thermal econcnfic model. He advocated the use cf zone 
analysis, as develcped by Wismer in Ref. 10, fcr large 
ccmplex systezs. In Wismer's approach to zcne analysis, the 
system is Ercken dcwn into zones containing one cr more 
ccmpecnents. The zcnesS are connected by linking variables. 
Twe methcds cf optimization are fossible. For the first 
methcd, in eé€ach zone, the zone parameters are optinized 
while the lirking variables are held fixed. fhen, the zone 
parameters are held fixed while the linking variakles are 
Optimized. Tkeé process iterates between optinizing the zone 
parameters and the linking variables until convergence 
criteria are satisfied. In the second methcd the zcnes are 
cut off frem each other and each linking variatle is 
separated intc two variables, one on each side of the 


cutting flare. Tke optimizaticn proceeds until the linking 
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variables cn e€ach side of the cuts are equal. 


Sheprard's system did not produce vapor in a boiler Eut, 
instead, maintaineé the working fluid in the liquid state in 
tke warm se€awater heat exchanger and vapcrized the werking 
fluid, amzcnia in his model, ty expanding it through a 
"plack kex" called a vaporizer. His optimization analysis 
was confined to the zone ccntaining the warm seawater 
exchanger, the feed pump, and the warm water circulating 
Euar. 


E. CEJECTIVES 


The present analysis impreves the previous model 
analysis ty replacing the liguid-to-liguid heat exchanger 
ané the vaperizer with a boiler. The boiler model is tc be 
cf sufficiert detail so that a realistic design is 
represented. The cther gajor components of the system, such 
as the turkire, the ccndenser, and the condenser circulating 
Fuaps, are included. A nonlinear prcegramming technigue is 
used to perfcrm the thermal-economic optimization of the 
ccmpleted system. The research reported here 1s intended to 
show the ctcsefulness of optimization analysis in the ¢esign 
of a lew tkermal difference system and te demonstrate some 
cf the prcblezs that must be considered in this tyre of 


analysis. 
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il.  Z2HE SODEL 


- => 2 Se =e 


A. GENESB2L LESCKIETION 


Vapcr peoweér generating cycles usually ccnsist of four 
processes. Fnergy is absorbed from a high temperature 
source and is used to vaporize a working £lugd ina “beiler. 
“Next, work is extracted by expanding the vaper in an engine, 
exhaustince tke vapcr at a lower temperature and pressure. 
The working fluid is then condensed back to the liquid state 
by rejectirg energy to the lcwer temperature energy sirk in 
the third recess. Finally, the liquid is returned to the 
higher presstre of the boiler by a pump. In the fcllicwing 
analysis, tke energy source is a current of ocean water 
Feated by tke sun and the energy sink is assumed to he 
cclder ocean water from deep currents. The scrking fltid is 
assumed tc ke propane although many other fluids are 


possible, such as ammcnia or the Freons. 


Six ccm@pcnents are modeled in the follcwing analysis. 
The hot side heat exchangers, the boilers, vapcrize the 
working fluid which is expanded through a turbine exhausting 
to a ccnderser that rejects the unused heat to the cold 
Sink. Three fumps are included, a feed pump returns the 
working fluid to the boiler and two pumps circulate the warn 
anc cold waters through the FEciler and the ccndenser, 
respectively. The arrangement of the system is shown ir Fig. 
2. The ncd¢s are numbered to correspond with the fluid 


state points shown later, in Fig. 3. 
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The analytical model consists of the equations 
descrikinc tke system and a set cf farameters to _ he 
cEtimized sc as tc minimize the capital cost of the flant 
subject tc certain ccnstraints. The objective function, the 
functicn tc ke minimized, is the sum of the capital costs of 
the pumps, the kEroiler and the condenser. Althcuch the 
turkine is a significant capital cost and would affect the 
design cf the system, it is not included in this model. 
There are twe classes of constraints. The first class is the 
explicit ccnstraints that specify the upper and lower fEcunds 
on the parameters. The other class of constraints is the 
implicit ccnstraints. These are the engineering equations 
that descrike the operaticn of the power plant. For 
example, a censtraint might say that the tegzperature inside 
the boiler must be less than the temperature of the warn 
seawater kecause the equations fail if the requirement is 
net satisfied. Ancther example would be a power balance 
that reguires that the sum of the energy flcews into and out 
of the syste ke zero. Additionally, they proscrike the 
regicn that the cptimizaticn rcutine may search within for 
the minigun cf the problem. The implicit ccnstraints are of 
two types, equality and inequality constraints. Inequality 
ccnstraintse set one-sided bounds on the feasible region. 
Scre of tke inequality constraints are restrictions flaced 
on the maximus OF tfinimum values of the farameters dte to 
judgements uade cn the basis of practical engineering 
consideraticrs. The equality ccnstraints are relationships 
that must Fe strictly satisfied at the optimimum. The system 
of eguaticns describing the model fit inte the fcllowing 
fcrm: 
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tirimize f£ (X) 


sumcet tO: 9 (4) 2 0, 1=1,2,.-4.2,0 
i 


h (X) QO, 1=m+1,....,0+n 
i 


where X is tke vector of the parameters tc be oftisized, 
£(X) is the capital cost function, g; (X) are the inequality 
ccnstraints and hj (X) are the equality ccnstraints. The 
parameters are the dimensions of the kEcilers and the 
concensers, the boiling and condensing pressures, the sea 
water velccities, and the mass flew rate cf the werking 
fivid thrceuch the systen. 


The kasis fer the economic framework pust be set cut at 
the start. Ail measurable costs are relevant to the final 
design cf tke power plant, but, for some types of analysis 
and design, many ccsts can be ignored. The life cycle cost 
of the CTEC must be considered before the decision is made 
tc commit the U.S. and the rest of the world tc the 
ccnsegquences cf chcosing this energy generation scheme as 
even a partiél soluticn to the energy prctlen. The life 
cycle ccst ccnsists of the capital costs, the financing 
costs, and the operating costs. The trade-cff of cperating 
ccst and capital cost is not considered in this resé€atrch, 
because the cperating cost 1s, for the most part, 
independent cf the overall dimensional characteristics of 
the compcnents. Operating costs are dcminated by the 
structure, site location, and maintenance requirements. 
Such things as the kind of instrumentation, ccntrols, 
kearings, and structural materials nave mcre effect cn the 
cperational ccsts than the length of the Eciler or the 
diameter cf the foiler shell. Therefore, cnly the capital 
cests of the majcr components, except the turbine, are 


considered. 
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To be tseful for model analysis cost data gust be 
transformed into a cost estimating relaticnship (CER) that 
shews hcw cost varies with some variable or some combination 
o£ variatles. The costs must be demonstrated to Le highly 
correlated with the variable chosen. Since the cest 
relation is cnly a correlation, the bcunds on the accuracy 
of the eguaticn shculd be specified by those who develce it. 
Since the relation is arrived at from a finite body cf data, 
the CER is directly applicable only in the range of the data 
ccllected. Eatrapolation cf the CER outside the rarge is 
dargerous recause there is nce data in that region tc support 
the assumed curve. However, the value of CER's is the use 
of them tc enakle the estimator to make a ketter analysis of 
a new déesicn than he could if the informaticn is not 
availartle. The relevant range of the CER must be keft in 
mind wher analyzing the solutions from the standpoint cf the 
reliability cf the answers. Another factcr tc consider when 
using cost estimating relaticnships is whether the equipment 
the data was taken from 1s ccmparable to the systen the 
analyst is considering. If the design, the technclogy, or 
the applicaticn of the proposed equipment is significantly 
different fircm that from which: the data is taken, tke CER 


may give unreliable answers. 


This aralysis begins with the develcpment of the 
encineerirg equaticns that ccnstrain the system. This is 
fcllowed by the development of a cost framework kbtased on 
capital costs as the function to be minimized. The 
engineering and économic relations are linked togettker in 
the cbhjective functicn and the constraints. A description 
of the algcrithm that performs the optimizaticn is included. 


B. CYCL& ANALYSIS 
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The tuilcing of the model begins with the analysis of 
the therocdynamic cycle. Figure 3 is a temperature-exntropy 
diagram cf tke werking fluid on which the cycle is shown. 
These nugrbered nedes represent the thermodynamic state 
points at that location in the system and ccrrespond tc the 
numbers in Fig. 2. The boiling and condensing process€és are 
assumed tc ke reyersible processes. In the pump anc the 
turbine the reversible process is indicated by the subscript 
tst' on the state pcint number. In the actual systenr there 
are irreversitle losses inside each component that are left 
out of tke analytical model or lumped in with cther Ilcsses. 
The pressure drops due to frictional losses in the tEoiler 


and the ccndenser are negiected. 


State pcint 1 assumes that the liguid leaving the trciler 
is saturatec. State points 2s and 2 represent the 
conditions leaving the feed fump. If there were no lcsses in 
the pump, the fluid would be at 2s. The fluid enters the 
beiler tc te heated, first, to saturaticn conditicns at 
state point : then boiled at essentially ccnstant fressure. 
State pcint 4 shows the vapor leaving the boiler as a 
Saturated varper. The lines from point 4 to 5s and 5 
represent tke expansion through the turbine. The path tc 
state feoint 5 acccunts for the losses due tc frictior and 
cther irreversikle processes in the turbine. The remaining 
heat is then rejected from pcint 5 to point 1, again assumed 


to ke a ccnstant pressure prccess. 


The tkecretical heat transfer rate in the boiler and the 
ccendenser is found by subtracting the enthalpy at state 
rcints 2 and 1 from the enthalpy of state points 4 ard 5, 
respectively, and multiplying by the total mass flow rate of 
the working fluid. 


Qab = Muse (Hy — Ha) (1) 
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FIGURE 3. THERMODYNAMIC CYCLE 
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Qrej = Mug ls - Ay) (2) 


The enthalpy at state point 5 is determined by first finding 
the moisture ccntent of the fluid exiting the turbine. 


Su- Ss, 
fie) = can ale 
Sta 


The enthalpy at point 5s is determined from 


Hc, = H, +2: Beg (4) 


This allcws B, to re found fron 


Hae 
ale anls {5) 


lee = 
Hope te 


where the efficiency of the turbine is assumed kKncwr. The 


power out of the turbine is now known. 


Peun = Swe(Hy ~- Hs) (6) 
For the feed pump the change in enthalify is accounted for by 


the change ir pressure assuming negligible rise ir the 
temperature. The enthalpy of point 2 1s determined frog 


UPR = Le opicmen 2b (7) 


where the efficiency is assumed known. The power required to 


drive the fugap is determined from 
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Pen = Myg (Hy - Hy) (3) 


fc gain a better feeling fcr the low thermal fetential 
systers end tc make comparisons with other systems several 
efficiencies are worthwhile calculating. The first is the 
Carnot efficiency which defines the maximum efficiency 
Fessible for a system operating between two temperature 
linits in decrees absolute, °R (°K). 


Teold 


lam = 1- 
Thot 


(9) 


The next efficiency calculation takes into acccunt the 
theoretical thermcdynamic cycle and is the ratio:of the net 
power cut cf the cycle to the heat taken up Ey the cycle 


Piur - P 
Ceye = oe (10) 
ab 


This acccunts for the thermodynamic inefficiencies inkerent 
in the cycle. Finally, the system efficiency takes into 
acccunt the cther internal consumers cf power: the 
circulaticn pumps, the pipe losses, and any other lcsses, 
such as tke antifouling equipment. This calculation 
acccunts Da(op g ail of the mechanical and electrical 


inefficiercies of the system. 


NET POWER AT THE EUS EAR 
Qab 





sys = (11) 


For example, if the hot water is at 75°F (23.9°C) and 
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the ccld water is assumed at 40°F (4.4°C) the Carnot 
efficiency is 6.5%. For one typical set of assumpticns, as 
shewn in Table 5, where the turbine and the pumegs are 
assumed te Fe €5% efficient, the cycle efficiency is 2.6% 


and the system efficiency is 2.4% for one feasible design. 


C. ENGINEERING FRAMEWORK 


This secticn describes the design of the components that 
make up the system. The level of detail of the analysis 
prcevides a physically realistic model and clcses the lcog. 
Tké pipe friction licsses on the cold and hct water inlets as 
well as tke fife lesses in the working fluid side are 
neglected. The turbine is included only as an energy 
transducer and no attempt is made to obtain costs, fhysical 
dizmensicns, cr fluid flow characteristics. fhe efficiencies 
of the turbine and the pumps are assumed as known ccnstants. 
Ike heat akscrbed in hneating the working fluid to the 
saturated state is a small pertion of the tctal energy taken 
wp ain tke kEciler; therefore, the sensible heating cf the 
fluid is neglected. 


1. TIke Eciler 


The Eciler ccnsists cf a number cf tube-in-shell 
heat exchangers connected in parallel as shcwn in Fic. 2. 
Figure 4 shows the internal arrangement of cne shell. The 
tute burdle is horizontal and boiling occurs on the outside 
surface cf tke tubes. The tube bundle is sukmerged ir the 
werking flui¢. The working fluid is maintained at a level to 
just cover the top rew of tubes. The warm seawater 1s fumped 


through tke inside of the tubes. 
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Ike clearance between the tubes is fixed at Q.5 
inches (1.27 cm) to allow for flaring of the tube ends or 
cther space and strength requirements. Similarly, the tube 
thickness is fixed at 0.028 inches (0.071 cm), since no 
strength reguirements due to differences kEetween internal 
and external fressures on the tubes or material loss due to 
ccrrosion of the tube surfaces are accounted for. 


At the tcp of the tube bank the leaving vapor 
ceontains a large amcunt of liguid dropiets which must be 
elizinated eitker in the boiler or by an external separator. 
In the mcdel the separation takes place in the bceiler by 
allowing sufficient space aktove the tube bank. Palen and 
Spall, Ref. 11, presented a methcd for finding the necessary 
sé€paraticn vclume that considers the effects of surface 
tensicn, density of the vapor and the liguid, and the mass 
ficew cf tke wcrking fluid. First the allowable vapor lcad is 
found frcs 


~~ | a Ce5 ; 
L = 22902 ae (12) 


Vagper loaé is an expression for the mass flow rate cf vapor 
threugh a urit volume. The expression for the maximum vapor 
lead limits the flow rate so that the moisture droplets are 
able to settle out of the vapor before it leaves the bciler. 
Ther the reguired segmental area is found frer 


4 
VL LS 


The segmental area is the cress sectioral area from the tor 
of the tuke bank to the top of the shell. Kncwing the 


segmental area the height of the tube tank above the main 


Se 





diameter is fcund fron 


2 2hs 
SA = €.25 [Ts ARCOS 
Ds 





Z 2025 
, bs (0.25eBS sls ) (14) 


Palen and Small also present a much simpler method of 
determining the total height of the tube bank. The kEeight 
of the separation vclume is set at 40% cf shell diameter 
which results in "hs" being equal to 0.1 Ds. Either trethod 
can be used alone or the first method can be used as a 
constraint on the seccnd. For this fproklem the segzrental 
area is fcund by bcth methods with the first used as a check 
on the values given Ey the second. 


Since the heat flux calculation is based on a single 
tuke, the nuxker of tubes per shell must be determined in 
order tc find the total energy aksorbed by the heat 
exckanger. The tubes are assumed to be laid out in a 60 
degree triangular pattern for maxiwum ccmpactness. The 
tukes are ilecated at the vertices of an equilateral 
triangle. Iwce cross sectional areas are calculated. The 
area of a Léxagon whose minor radius is equal to cne half 
the pitck is determined fron 


P\2 
A = «(5 TAN 30° (15) 


Then the effective cross sectional area of the tube bank is 
fcound by subtracting the clearance area, SA, from the cross 
sectional area of the shell after correcting the shell 
diareter fer a clearance space to aliow for free flow 
channels both within the tuke bundle and arcund the ottside 
of the tute Eundle. The effective cross sectional area of 
the tube cLark is divided by the area of the hexagen, found 


ES, 





akcve, tc determine the number of tubes in a shell. 


The flcw of the liguid and vapor through the tube 
bank is assured to be turbulent. The inccming subcccled 
working fluid fren the feed pump sixes with the 
recirculating fluid returning frem higher uf in the tube 
bank. fhe circulation of the fluid is induced ky the 
boiling process. This mixture enters the tube cank and 
begins the krciling process. Near the bcttom of the tube 
bark the Feat transfer goes mostly to Ering the fluid up to 
the boiling temperature. After the fluid reaches satturated 
liguid ccnditicns, stable boiling begins. The vapor bubbles 
rising threugh the surrounding liquid increas¢ the 
turbulence arcund the tubes. As the liguid continues uf 
thrcugh the tube bank, the local boiling temperature is 
decreasing due to the decreasing hydrostatic head. This 
decrease in the boiling temperature increases the 
temperature difference and increases the heat ficw and 
increases the bubble generaticn rate. The vapor generation 
may recome sc rapid that a condition called vapor blanketing 
may ccctr in the interior of the tube bank where the liguid 
cannot flow inward fast enough to displace the vapor flowing 
through the region. This is similar to opurnocut ina 
cenvective Eciling in that the tubes dry out and the heat 
transfer rate decreases. This effect becomes proncunced 
cniy near the critical heat flux according to Starczewski, 
in Ref. €. 


All cf the heat and grass flow variaticns in the tube 
bundle are averaged by performing the calculaticns for a 
Sirgle tuke lccated at the mid height of the tuke bank. 
This tute is assumed to represent the average of all the 
tukes in the tube rank. The boiling situaticn is assured to 
be saturateé fool toiling. Quoting frcem Ccllier, Ref 12, 
"rool poiling is defined as boiling fren a heated surface 


sukmerged in a large volume of stagnant liguid". Since the 


34 





fluid vclume seen by the tube is neither large nor stacnant, 
a ccrrecticn is needed for the effect cf being in 2 tube 
bundle. Saturated boiling is defined as boiling at a 
constant temperature and pressure as the liguid changes to a 
vapor. This is frobabiy not strictly true because the 
boiling temperature is ccntrolled by the lcecal presstre at 
€ach tuke, and the fluid is not reaching a stable 
temperature as it rises in the bundle. The situation is 
prckakly closer to convective boiling since the fluid is 
boiling in a confined channel between the tubes. An 
additioral effect in the kEkoiler is the fact that the 
hydrostatic read at the bottom of the tuke Eundle causes the 
boiling te take place at a higher fressure at the tottom 
than the kciting at the top of the bundle where the pressure 
is equal tc the koiler exit pressure. This effect 1s 
acccunted fcr Ly using the average head of the bandle in the 
heat transfer equaticns. 


The heat transfer analysis of the bciler begins with 
tke kasic eqvation 


Q = UA AT, (16) 


where U is the overall heat transfer coefficient, A is the 
area thrcugkt which the heat is transferred and AT, is a 
prceperly defined mean temperature difference. In, Sens 
paper, the prceduct UA is treated as a single guantity called 
tke thernual ccnductance. One may look at the equaticn as an 
analogue cf the basis electrical eguation I = E/R where I is 
equivalent tc Q, E is the same as AT,, and UA is 1/8. This 


analogy is used when developing the expressicn for JA. 
In crder to account for the effects of the tube 


seeing many nearby neighbors, the equation is modified by 
tke incltsicn of a factor, F, to account for the effects of 
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the surrcundirg tukes. Equation 16 is changed to 


Q = UA F AT, (17) 


where F is defined as the ratio of the heat flux of the tube 
Eurdle tc the heat flux for a single tube in a stagnant fool 
boiling situation. According to Ref. 13, a report published 
ky werkers at Heat Transfer Research, Inc. (HTKI) of 
Alhambra, California, F is greater than cne. In Figure 5, 
which is sketched frem a figure in their report, F franges 
from 4.72 at AT,, = 8°F to 2.63 at AT, = 40°F. Since ne data 
is given cn the conditions of the test, F is assumed to he 
fixed at three. 


Ir crder to tind ATy. the heat exchanger 
effectiveress, €, is determined. The effectiveness is 
defined as tke ratic of the actual heat transfer rate to the 


Maxipum fcossible heat transfer rate as shown by 


Ch (Thin ~ Trout) 


i aaenEEEEP Eee (18) 
Cmin (yin ~ Tein ) 
where 
C= (M Cp) tiuid (19) 
Cerin = Minimum of Cy or Cy, (20) 


In the bciler, Ty jin is the temperature of the entering 
Seawater, ard Tein is the saturaticn temperature cf the 
working fluid at the pressure at the mid height cf the 
bkeiler. Reference 14, by Kays and ILIcndon, cortains 
eguations fcr finding the effectiveness fcr varicus flow 


configuraticrs. In addition to finding Cai, Cea muse ee 
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FIGURE 5. COMPARISON OF A SINGLE TUBE 
VERSUS A TUBE IN A BUNDLE 
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kncwn in crdéer to decide which equation afrlies. For the 
case of Eciling, Cy,, is infinite; therefore, regardless of 
the ficw pattern 


€ = 1 - EXP(-NTU) (21) 


where NTC, the number of heat transfer units cf 2a heat 
exchanger, is defined as 


NTO = UA/Cmin (22) 


and 


Cmin = Mew Cpsw (23) 


AT,,is found fren 


UA AT 
Cmin ATo 


and the temperature of the sea water at the outlet cf the 
tuke is detéerained frcem 


Cos THin ~ Trout (25) 
AT, 
where 
AT, = Thin7 fein (26) 
The quantity UA is determined next. Figtre 6 
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illustrates the sgethod used. If UA 1s pictured as the 
reciprcecal cf resistance then the flow cf energy froma one 
fluid tc tke cther can be pictured as a series electrical 
circuit and UA is determined fron 


DA = 1/(R1 + R2 + R3 + RY) (27) 


Reference 15, by Holman, shows how each resistance can be 
found. &k1ti1is the thermal resistance due tc convecticn on 
the inside cf the tubes, the sea water side. 


, 
Me = See (28) 
TT Dey dj; L 


R2 is the thergal resistance of the tube wall. 


7 in(d, /4d,; ) 


= aeeetiaaleee (29) 
27k, L 


R3 is tke thermal resistance attributed te boiling cn the 
cutside cf tke tukeés. 


R3 (30) 


min al, 


R4 is the thersal resistance due to any other resistances, 


such as fculing and ccrrosion defosits. 


{ke keat transfer coefficient on the seawater side, 
hye is determined from the dimensionless group called the 
Nusselt pumker defined as follows for flow inside of tubes. 
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Nu = ——— (31) 


In crder tc tind the Nusselt number, three distinct regions 
of flow are allowed for in the precklenm: laginar, 
transiticr, and turbulent flow. The divisicn points are in 
terms of valves of the Reynolds number, Re = udp/y, 


laninar: 0 < Re 2,060 
transition: 2,000 < Re < 10,000 
turbulent: 10,000 < Re 


lA 


lA 


ent ini oy 


In the lazinar region the Sieder-Tate correlation is used. 


173 | 
Nu, = 1.86 (Re Pr dj /L ) (32) 


In the ttrkulent region the Littus-Boelter correlaticn is 
used. 


or: ae 
Nuy, = 0.023 Re Pr (33) 


In the transition region the flew is unstable, but frog data 
in Ref. 14 it appears that a fairly smooth curve drawn from 
the lamirar curve to the turbulent curve adaguately 
represents tke value cf the Nusselt number in the transition 


regicn. Tbis is done by assuming the forn 


Nu-. = y Nu, + (1 - y) Nur, (34) 


where 


ee 


S 2 
y = aRe + ©— Re + Cc RE + @ (35) 
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The coefficients in the polynomial are found kL-y requiring 


Nuy.. =U eat eS 2,000 
Nut. = Nuy, at Re = 10,000 


and that tke first derivatives with respect to the Reynolds 
nugkrer match at Re egual to 2,000 and 10,000. 


The I.1L. Mostinski correlation, from Ref. 16, for 
single hcrizcrtal tube pool boiling is used to find the heat 
transfer coefficient on the outside of the tube. ATRI's 
extensive experimentation led Palen, Yarden, and Takcrek, in 
Ref. 13, tc conclude that Mostinski's equaticn produced the 


“acst consistent results". His correlation is 


C.69 0.7 0.17 Te 1C 
Bye = 0.00658 po 86g) (1.8 x +4r +1027 ) (36) 


where 
I = E/Pe 
Pe critical pressure of the working fluid (psia) 
g heat flux (BTUyhr-ft?) 


A final heat flux relation needed is the relation fer the 
Maximum heat flux allowed. Mostinski,s ccrrelation fcr a 
sirgle ture is used. 


0.35 0.9 
= 803 PR © (1 - x) (37) 


Imax 


where the units are the same as in eguaticn (36). From 
HTEI's data, under unknown conditions, gy., for the kundle 


appears tc ke cne third of the q,,, for a single tute. 
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Ike fewer reguired to pump the salt water ttrough 
the tubes depends cn the pressure drop across the tuke bank. 
There are three causeS of the pressure drop. The flow 
experiences a loss as the water enters into and exits fron 
the tubes in addition to the frictional loss along the tube. 
From Ref. 14, the pressure drep fcr a liguid is determined 
fren 


2 4£f£L 
Ee = 105) Ceuta | soe + Ko - Ke) (38) 
} 


In crder tc determine f, Ke and Kc, tke flcw regime must be 
specified. Tke regimes are a little different than those for 
the Nusselt rumber. . 


laminar: QO < Re < 2,0C0 
transition: 2,000 < Re < 5,000 
turbulent: 5,000 < Re € infinity 


The Fanning fricticn factor, f, in the laminar regime is 


f = 16/Re (39) 
In the turktulent regime, the Fanning fricticn factcr is 


found frcem the Blasius equation in the range 5,00C < Re < 
1€€C,000 


0.25 
f = 0.079/Re (40) 


which is a close approximaticn to the Karman - Niktradse 


egtaticn, 
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“90.5 G25 
(4 f) = -0.8 + 0.87 In(Re(4 £) ) (41) 


which is valid thrcughout the turbulent regicn (Re > 5,000), 
assuming smccth walled tubes. Either equaticn (40) or (41) 
may ke used with similar results but the Elasius is Letter 
Suited tc this analysis because the Karman-Nikuradse 
equation must be sclved hy iteration. Equaticn (41) 
consumes more computer time to solve and can cause accuracy 
problems if not enough significant digits are evaluated. 
The fricticr factcr in the transition region is determined 
by fitting ¢ least squares curve through a hand drawn curve 
frcem the end cf the laminar curve to the keginning of the 
turkulent cturve. 


The entrance effect, Kc, and the exit effect, Ke, 
depend on the Reynclds number, the free flow to frontal area 
ratio, $, and, in the laminar region, on the tube length. 
For the tEciler, $ is determined by 


= aS Se ee 42 
® ; (42) 


Figure ¢, fage 7-13, in Ref. 17 is used tc find Ke and xc. 
The effect cf Reynolds number and L/d ratio was eliminated 
from the relationships in the laminar and the turrulent 
regimes cf tke flow, tut the Reynclids number was included in 
the trarsiticn region because of the large variaticn of Ke 
with Reynclds numker. The fcllowing equaticns are used to 


find the entrance effect in the various flow regimes. 


Laminar: (using the curve for 4 (L/d)Re = 0.1) 


Ke, = 1. - 0.4 ty) (43) 


4d 





MMurbulent: (using the curve for Re = 5,00C) 


Key, = 0.52 - 0.4 > (44) 
igansiticr: 
Key, = A Ke, + (1 -A) Key, (45) 
where 
XN =e Re+ef (46) 


The ccefficients e and f are determined Ly requirirg the 
transiticn Ke to equal the laminar Ke at Re equal to 2,000 
and the turkulent Ke at Re egual to 5,000. The exit effect 
is the same fer all flow regimes, and is found fren the 
fcllcewing equation that approximates the 4(L/d)Re =0.05 
curve in Kef. 17. 


2 3 
Ke = 1.€0093 - 2.5178 + 1.1613 - 0.176776 (47) 


2. dhe Ccndenser 


Tke éesign of the condenser is almost the same as 
the desian cf the boiler. It is a fixed tuke sheet, 
tuke-in-stell heat exchanger. The condenser does nct reed a 
separaticn vclume so the entire shell is filled with tukes. 
In crder to allow the entering vapor to reach all cf the 
tutes, ckannels through the tuke Eank are needed. The 
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allewance for this 1s accomplished by decreasing the shell 
diameter by a fixed factor when calculating the numter of 
tukes. The calculaticn of the working fluid side heat 
transfer coefficient, h , is different. The equation 
acccunts for the ccndéensate from one tube falling cn those 
belcw it ard for the subcooling of the ccndensate. 4M.M. 
Chen's eguaticn, as developed in Ref. 18 and shown in a more 
Single fcrm in Ref. 17, is used. 
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CoAT = k H! 
Be = Co7ZE +0. ( P )ca=1) | [2A A, “Ay)% Bia) (48) 
Hf, pa de ff AT 
where 
k = BIU/ft2—-hr-lOF 
H' =H + : AT 50 
fg "tg * 9 © a 
n = number of tubes in a column 
HY is the entkalpy corrected for subccoling of the 


cerdensate. Egquaticn (48) replaces equaticn (36) ir the 
analysis cf the condenser. In the circular shaped tube 
bark, n is fcurd by determining the number cf tubes in a 


cclumn of average height. 
Ds 
n= 9.425 — - 1 (51) 
pa 


The definiticn of AT, becomes 
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AT, = Isat ~- Tswin (52) 


in contrast tc eguaticn (26) above. The seawater cutlet 
temperature is faund from 


a Tswout ~ Tsw,in 


53 
oe (53) 


which replaces equation (25). 
In crder tc find the heat transfer ccefrficient, the 
temperature cf the outside tube surface is needed in 


equation (49); therefore, the thermal netwerk is solved for 
the outside wall temperature and for the heat flow rate. 


Q = Ustt — ec) (54) 


where 


UA' = 1/(R1 + R2 + R4) (55) 


Toulk = average of the inlet and 
the outlet seawater temp. 


There is ne factor, F, in the heat transfer equaticn for 


ccordensation since the effect of neighboring tubes is 


acccunted fer in eguation (47). 


Trere are three pumps in the model. Poth 
Circulaticn fumps for the bceiler and the condenser are 
assumed to ke propeller pumps (axial flow pumps). The feed 
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FuRES are assumed to be centrifugal pumps. The differences 
in the characteristics of the two types of pump show uf cnly 
in the calculaticn of capital costs of the fpumps ir the 
model. 


Ike chcice of the type of pump to aprly is made on 
the basis cf the flow rate, GEM, and on the head the pump is 
working against. Cualitatively, the propeller pump is used 
where tke flew rates are high and the head is lcw. The 
centrifugal fumps apply where the flow rate is comparatively 
low and the head pressure is high. The qualitative picture 
is placed in tetter perspective through the calculaticn of 
the specific speed fron 


0 
ae RPM GEM _ 
ee 0.75 

Head 


The centrifuqal pump is used when 500 < Nsgq < 7,500 and the 
propeller pugp applies when 7,500 < Nsg $ 15,000, acccrding 
to Ref. 16. The choice of the type of pump is checked after 
the optisizaticn process has been completed by sclving 
eguation (55). If the value of the specific speed is 
outside cf tke range, 500 < Nsq < 15,000, the maxinun 
efficiency cf the pump drops off rapidly; therefcre, few 
fuarps are designed tc operate outside cf this Lange. 


Since the pumps are "Elack boxes", only the power 
consumed ky the pumps is calculated. The pewer required to 


drive the punop is found fron 


M Ap 


P= (57) 





Ppump 
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CL. COST FEAMEWORK 


As discussed at the start of this repcert, the capital 
costs of tke major components are the cnly costs to be 
considered ir this mcdel. The bciler and the condenser are 
both tuke-in-shell heat exchangers assumed to te cf the 
fixed tube skeet type. The only difference keing that the 
Shell of tke koiler is fartly filled with tubes. Both 
seawater circulaticn fumps are assumed to be propeller type 
Fugps and sade cf the same material. The size and number 
per heat exchanger shell may ke different fcr the bciler and 
the ccndeéenser, The number of working fluid feed fumps, 
which are the centrifugal type pumps, are assumed tc equal 
tke number cf Ecitler shells. 


The ccst estimating relationshifs iGE the heat 
exchangers ard the centrifugal fumps are taken from Ref. 20, 
by K.M. Guthrie. The eguation fer the cost of the heat 
exchangers is 


0.63 


Sep = 118-3 (Fp Ba), Pee teeta) (58) 


where the ccrstants are calculated from a Ilcg-log flict of 
cost versus tctal tube surface area and 


= 0.€ for a fixed tube sheet 
co 0.0 for pressure correction (p < 150 psi) 


raj 
a. 
' 


Be = patéerial factor from Table 1 
Ine = ccost index 
= dg LN 


fer a relevart ranse cf 100 < A < 10,000 fté. 


The ccst equaticn for the centrifugal pumps and the 
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asseccliated mctcrs is 


: 0.721 

$e, = 1.013 F, Fy Ig, Cy (59) 
where the constants are fcund from a straight line 
approximaticn cf the log-log plot cf ccest versus C for the 
section cf tke curve from C = 30,000 to Cc = 300,000 fer an 


electric motcr driven pump and 


Fo = 1.0 fcr suction temperatures < 250°F (121°C) 
Ig, = ccst index 

Fy, = Baterial factor from Table 2 

Cy = predtct of the flow rate in GPM 


and the pressure differential in fsi 


for a relevant range of 30,000 < Cy < 300,000. The cost 
Gata in Guthrie's article is given fcr a time kase of 
b1d-1968. | 


Reference 21, also by Guthrie, gives cost data for 
propeller pumps and centrifugal pumps. The equation fcr the 
prepeller fugofs with the motors is 


0.783 
Sep = 5-767 Fy Fy Igp Flow (6 0) 


where the ccrstants are found from the log-lcg plot cf cost 
versus flcw rate for an electric motor driven pump and 


Fp 1.0 for a suction pressure < 150 fsi 
Fim 
Iep = cost index 

Flow = flew rate in GPM 


material factor from Table 3 


for a relevant range of 1,000 < Flow < 100,000 GEM. The 


tine base fcr the frepeller pump data is fcr the erd of 
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Table 1 
BATERIAL FACTORS FOR THE HEAT EXCHANGERS 





Material 
Shell Lube Fo 
Carkcn Steel Carbcn Steel 1.00 
Carkcn Steel Brass 15.2 
Carkcn Steel Stainless Steel 365.2 
Stainless Steel Stainless Steel 4.50 
CarkEcn Steel Monel 3 
Ncnel Monel 4.95 
Carben Steel Titaniun Vi310 
Titariun Titanium 16.60 
Table 2 


MATEFIAL FACTORS FOR THE CENTRIFUGAL FEED FUFFS 


Material Fn 
Cast Iron 1.00 
Bronze 1.28 
Cast Steel Ve32 
Stainless Steel (Lass 
Monel Jeo3 
Hastelloy C 2.89 
Titaniun 8298 


Deu 





Table 3 
MATERIAL FACTORS FOR PROPELLER PUMPS 


pumps 
Material Fn 
Cast Iron 1.00 
Cast Steel 1.28 


Stainless Steel 1.64 
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1970. Ike costs arrived at from equations (58), (59), and 
(60) are in dcllars per component, for example, dollars per 
heat exchanger shell. The purpose of the cost incex is 
allcw all the equipment to be costed on a ccmmon basis. The 
tige difference is assumed to be small encugh so that the 
cost index is set to 1.00 for all the components. 


E. ASSEMELY CF THE MATHEMATICAL MODEL 


The engineerirg and the ccst relationships are linked 
together in the model. The cost equation is the function to 
be minimized and the implicit and explicit ccnstraints are 
the engireerirg equations that restrict the feasible region 
of the ercetlem. The explicit constraints restrict the 
Variables, tke vector X, to positive values. The rest of 


the model is fortulated in the fcllowing fcrs 


minimize £ (X) 


subject to: g (X) 2 0, i=1,2,....,a 
a 


0, L=m+1,....,mMtn 


h. (x) 
a 


where tke definiticns of the elements of xX are shown in 
Takle 4. Several of the variables must have only integer 
values if the system were actually built, but for the 
objectives cf this aralysis they are left as continucus real 
variables. Fer example, the number of shells may have an 
integer value in the actual system but this requiremert can 
ke accomccated by solving the system for a real value then 
integerizing the number of shells and solve the prckleg with 
the next higker and the next lcwer integer number of shells. 
The lower cf the two solutions would be chosen. The final 
rescluticn cf the correct number oF shells would 
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x (1) 
X (2) 
x (2) 
x (4) 
x (5) 
x (6) 
x (7) 
x (&) 
x (9) 
x (1C) 
x (11) 
x (12) 
x (135) 


it 


Table 4 
DEFINITION OF THE ELEMENTS OF THE VECTOR X 


velocity of seawater in the bciler 


cutlet pressure 
inlet WY 


it " condenser 
in the boiler 


it " condenser 


tuke length in the boiler 


mass flow rate 
diameter of the 


“ tI Te 


outside dia. of 


" condenser 

of the werking fluid 
boiler shell 
condenser shell 

a tube in the boiler 


Lt u " condenser 


number of bciler shells 


at " condenser shells 
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prckably defend on other than thermodynamic and cost 
consideraticns, such as a requirement that cne extra shell 
be included to allow for one shell teing down for 
Maintenance. Chviously, the number of tubes in a shell must 
be an irteger value but one extra tube cut of several 
thousand is a ginor matter. In additicn, there may be a 
reguirement fcr extra tubes to allow for a percentage tc be 
plugged wken they develop leaks so that the entire shell is 


net put ctt cf ccmfission by a minor leak in one tute. 


The ckjective functicn is the sum of the costs of 
the boiler ard the condenser, and the punmgs. The model 
assumes tbat there is one circulating pump and one feed fump 
for each stell in the kciler and that there is one 
circulatirg fump for each shell in the condenser. Since 
there is ne cost analysis for the turbine, nc assumfticn is 


needed akcut the number of turbines in the systen. 


There are Nineteen constraints in the present 
prcgram, cne cf which is optional. The first three ccncern 
the saturaticn pressures in the boiler and the ccndenser. 
The next twelve place restrictions on various dimensicns of 
the Eciler and the condenser. The last three require that 


scre intermediate calculations stay positive. 


Csing the cbhjective function as a constraint has 
proved advantacecus in cases where the probleg was 
ncen-ccnvex or where the program had difficulty in finding a 
solution. It is used by setting a boundary cn the chjective 
function <c that the constraint is either infeasible ocr 
Slightly feasible at the starting point. The proper chcice 
must be icund by trial and error. One caution to note is 


when equality constraints are involved. If the starting 
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Fcint is near the optimum, the program tries to move away 
frcm the bcundary and the equality constraints may not be 
satisfied wken the pregram finishes. 


Tke need £or the cther eighteen inequality 
ccrstraints is net obvious when the preblem is first 
formulated. The thermodynamics of the proklem require, for 
example, that the fropane saturation temperature Fe less 
than séatater temperature in crder for the Eciler te aksorb 
energy, crt the pregram has no way of knowing that, unless 
it it specified. If the program can imprceve the cost and 
satisfy tke cther constraints Ey reversing the relaticnship 
it will dc so. This inversion may take flace only fora 
shert tige during the sclution of the problem, but, if the 
equations fail under this ccndition, the ccnstraint cn the 
temperatures in the bciler is necessary. Mest of the other 


constraints are needed for the same reason. 


Scme ccnstraints are net added to the probler since 
they are net binding in the solution. There are twe such 
ccnstraints. The first is the segmental aréa required to 
allow sefaraticn of the liquid droplets from the vafer. The 
segmental area set aside by the 40% rule-cf-thump is more 
thar that reguired by the consideration of the physical 
preoperties cf fprepane. The cther ccnstraint is the 
restrictica cn the maximum heat flux by Mostinski's 
correlaticn, eguation (37). The constraints do not Eecome 


binding in the scluticn. 


The three equality censtraints contain the 
calculaticns cf the properties of the various compecnents 
modeled in this analysis. The first cnée concerns the 
bciler. It states that the actual energy akscrption rate of 
the Eciler stells gust be equal te the energy aksorption 
rate creguired by the thermodynamic cycle. Similarly, the 
second eguality ccnstraint states that the heat rejected by 
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the condenser must equal that required by the thermodynamic 
cycle for the ccndensation of the vapor from state point 5 
back to state pcint 1 as shown in Pig. 3. The third 
constraint calculates the pcwer required by the various 
Fuaps and the pcwer extracted from the turbine to run the 
system's furps and to produce electrical power for 
consumpticn. This ccnstraint requires that the power cut of 
the turkine equal the power used by the fumps and the 
generator. The net useful power cut of the syste is a 
design cksective and is specified as 25HW herein. 


Since the arrival of the electronic computer, a 
great zany nonlinear programming techniques have teen 
developed frem early theories and more current treseéarch. 
The ccmputer program, called the Sequential Unconstrained 
Minisizaticn Technigue (SUMT), developed by Mylander, 
Hcolens, and McCormick in Ref. 22 is used to optimize the 
mcdel. Their prcgram implemented much cf the theory 
contained in Ref. 23, by Fiacco and McCcrmick, about 
ncnlinear programf#zing using unconstrained minigzization 
technigues. SUMT-Version 4 is chcsen because of its akility 
to solve a wide variety or froblems. The thecretical 
requirements for SUMT to ke guaranteed to find a local 
Minimum is that there exists scome point which satisfies the 
ineguality ccnstraints and that there not Le a local sininun 
at points where X goes to infinity. If certain convexity 
conditions are satisfied, then the local rinimum found is a 


glcbhal rinimun. 


SUMI uses a penalty functicn tc transferr the 


constraired cptimization proklen 
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finimize f (X) 


sUrmect tc: og) (x) = 0, 1=1,2,. 5-7 
z 


h (X) = 0, i=m+l,....,m¢n 
1 


intc a seguence of unccnstrained problems cf the forn 


m m+n 
mieye) = £{X) - © bln g. (X) + 17e..k -(h (y= (61) 
1=1 1 l=m¢+1 2 


where F(X,r) becomes the function to be sinimized. The 
inequality ccnstraints, g(X), cause a large fenality to cccur 
as a bouncary is appreached. The equality ccnstraints adda 
large fperalty whenever h(X) departs from zero in €ither 
directicn. The size of the penalties due te g(X) and _  4h(X) 
is contrclled by the arbitrary parameter "rr", The algcritho 
solves the preblem by solving a series of sukproblems. For 
Ty 20, >03--->0, SUMT minimizes P(X(r,),r,) in each sukproklen 
in crder to find the minimum of £(X) as "xr" tends tc zero. 
When a minimum of a subproblem is found, "r" is reduced by a 
factor assigned by the user and the new subprcblem is 
sclved. The Frocess continues until the convergence 
Criteria chcesen fry the user is Satisfied. Figure 7 
illustrates schematically the results cf the segquerce of 
sukprcbleas and the qualitative shape cf the . fenalty 
function. Three methcds of mpinimization are included in the 
Frcgtranm. Tke user may choose either the generalized 
Newton-Raphscn, steepest descent, or McCorpick nodification 
of the PFletcher-Powell method. 


There are several other options that contrel the 
operation of SUMT. The user may chccse the initial 
penalty, r, cr, if the problem satisfies certain conditions, 
the user may elect to have SUMT compute the initial "r", 
There are three choices for the convergence critericn cf the 
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P(X,rp) 





FIGURE 7. PROBLEM SOLUTION SEQUENCE 


Dh, 


INFEASIBLE REGION 





prcblem. The criterion used depends on the ineguality 
penalty terms keing close enough to zero. Ancther oftion is 
the choice of the calculation of the sukproblem stcpring 
criterion . Two of the choices compute the inverse cf the 
seccnd partial matrix and the third choice ccmputes cnly the 
gradient of the penalty function ( V,P(X,n)). One cf ethic 
mcst useful cfptions is extrapolation. After two or three 
sukproblegrs have been solved SUMT is able tc make a first or 
seccnd crcer extrapolation to the starting cint for the 
next sufprcklem. This can speed up the sclution cf the 
proklem meéeasvureably. 


Ike characteristics of the proklem limit the choices 
for each cpticn. The penalty factor couldn't ke set Ly SUMT 
due to the presence of equality constraints and because the 
initial ¥ is usually not close to any boundary. The stcprping 
criterion fcr the problem is that the ineguality fenalty 
term beccme less than some value, usually VOR. 
Extrapolaticn is used, but it makes little difference in the 
answer cr tc the sclution time whether a twc point or three 
pcint estimate is made. The subproklem convergence 
criterion is the magnitude of the gradient of £E(X,r) 


bkeccming less than some small number. 


Tke choice cf starting penalty, Ey affects 
significartly the solution time and even whether a solution 
is found ir scme frobklems according to kef. 7. If the 
initial ‘tr" is toc small, the subproblem crtimum may fre the 
pretlem's cptimum which may be hard to find cr require more 
caiculaticns than if a much larger "r" is chcsen. If "r" is 
toc large, tbe first few optima lie near the center of the 
feasible region and are unrelated to the optimum cf the 


oksective functicn. 


Scaling of the variables, constraints, and the 


cbhjective furcticon is necessary. The variatles must be 
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within two cr three orders of magnitude of each other. For 
example, the working fluid flow rate is of the order of 107 
whereas the tube diameters are of the order 109. If the 
finite difference interval is 10-*, the partial of the 
eguations with resfect to the mass flow rate may be lost due 
tc the rcundcff error providing the only Significant 
figures. With the constraints and the objective function, 
unscalin¢g may result in one term in the fenalty furction 
being dcminant, and the algorithm finds that most cf the 
decrease in the penalty function takes place in the dominate 
tern. The result is that SUMT works hard to decrease that 
term and does not see the other terms resulting in a 


meaningless sclution. 


In making first and second order estimates of the 
scluticn, SUMT may generate negative numbers in the X 
vector. When it evaluates the constraints and objective 
functicn for -this point, the eguations fail and may cause 
the program tc fail depending on the severity of the failure 
and how many times the error has cccured. This prcblem is 
circumvented ty testing the vector upon entry inte the 
sukroutine where the objective function and the constraints 
are evaluated and setting the ccnstraint Eeing evaluated 


equal to scme negative constant. 


The main program of SUMT is altered so that the 
prc-lem can be restarted at the current soluticn without 
reading in a new problem. This 1S necessary because the 
prcgram would give all the indications of having fcund a 
sclution, tut if the problem is restared it cculd often find 
a lew valve for the cbhjective function. The new soluticn is 
ccmpared to the previous solution and if the two solutions 
are net clcse enough the problem is restarted after making 
scme changes to particular parameters, opticns, Or 


ccnstraints depending on the objective the user has in mind. 
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LIL. RESUDPS AND CONGLUSIGHS 


282 GG GE a aS eS = SS a See ee 


One cf the majcr tasks in the development of a nonlinear 
mcdel is the linking of the model with the optimization 
prcgram in crder te analyze the model. The selecticn cf the 
prcper ncnlinear programming technique is essential tc the 
seclution. Several techniques require that the ckjective 
function and the ccnstraints be in a particular fcrm. For 
example, Gecmetric Programming developed by Cuffin, Peterson 
and Zener, in Ref. 24, requires that the equaticns be 
frosyncmials which is not true of the model developed ir this 
paper. SUMIT was chosen because cf its wide applicaticns to 
proklems cf any tyre. It only requires that there exist 
some feasible point and that local minima not exist at 
Per inity. It has successfully solved Dany ncenconvex 
proklems. 


Two versicns of the analytical model are presented. The 
first is called the iterative versicn because the 
calculaticn cf the heat absorpticn rate and heat rejection 
rate require iterating on some quantity to find then. 
Mostinski's correlation, equation (36), contains the heat 
flux which is an unknown. It is not possible to solve the 
system cf equations explicitly for the heat flux; therefcre, 
a SUCCESSIVE approximaticn technique is used. The 
calculaticn cf the heat rejected in the condenser involves 
the same process with the temperature of the outside tube 
wall as the iterating variable. SUMT is urabtle to arrive at 
a ccnsistent sclution in the time allowed, up to 90 minutes 
of computer time. The soluticns are quite different for 
each try ‘vith nc pattern.evident akcut tke directicn the 


ecluticn right lie. 
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There are several possible explanaticns fcr this 
behavior. FEerhaps the behavior of this version cf the 
prcblem is such that SUMT is not the best choice fcr the 
optimizine algorithm or possibly the prceklem cannct be 
sclved at all in this form. Another source of pfrecklems 
could be the numerical differentiation procedure. Numerical 
differentiaticn is subject to severe truncation and 
reund-off errecr. The subproblem always terminates when SUMT 
is unable to reduce the penalty functicn during the next 
iteration, frcm one estimate of xX to the next estimate. 
Remember the ccnvergence critericn, never satisfied for this 
versicn cf the problem, is the magnitude cf Vy E (XT) 
beccming less than some small number. No ccenclusive results 
could te cktained with this version so the second version is 
used fcr the present analysis. 


In this second version the heat transfer coefficierts on 
the FEropene side of the tube are fixed. This enakles the 
€guality ccnstraints for the bciler and the condenser tc he 
calculated to the full accuracy of the computer without 
increasing the time required. The ccrrelaticns, equation 
(36) and (48), are used for ccmparison with the fixed value. 
If the calculated value and the fixed value are within 20 to 
25% of cne another, the solution would be considered to be 
gccd sirce this is the tolerance claimed EOE the 


GCerrelaticns. 


SUMT seemed to be able to solve the new prceblem kecause 
a series cf seven soluticns have been preduced where the 
values of the elements of X and cf the objective furction 
vary nce more than 3.7% and 1.2%, respectively, akout the 
mean when different starting points and different SUNT 
contrcl farameters are selected. See Table 6 for a Summary 
of the scluticns ottained from the assumed data, as shcwn in 
Takle 5. The design implications of these results are quite 


reascnaktle. The range shows the spread in the scluticns. 


63 





Table 5 
PROBLEM INPUT 


Fouling Thermal Resistance 
Eciler and Condenser 
Elevaticn cf Bciler 
AEcve Condenser 
Channel Allowance 


Efficiencies for the Turbine 
and the Pumps 

seawater Temperatures - Hot 

- Cold 

Net Output Pcwer 

Eundle Factor (F) 

Ccest Irdexes (I) 

Tuke Wall Thickness 

Cistance Between Tubes 


Min. 


Heat Transfer Coefficients 
fcr Eciler and Condenser 


Ecoiler and Ccndenser 


Feed Furs 
Circulation Fumes 


64, 


DATA 


0.0 hr-°F/BTU 
(0.0 °C kh) 
158 Et 

(4.57 m) 

0283 £¢ 

(0.25 m) 

85% 


or 

(23 .9°C) 

4OQOF 

(4.4°9C) 

8.532x10? BTUyhr 
(25 MW) 

3.9 

1.00 fcr all ccmp. 
Q.028 in 

(0.071 cn) 

0.5 in 

(1.27 cm) 

145 BIUyfte-hr- OF 
(0.822 kW/m2-°C) 


Carbon Steel Shells 
Brass Tuces 
Bronze 


Stainless Steel 





ariacle 


Bciler: 
Velocity 


Eressure 


fluke length 


Shell dia. 


Tube dia. 


# of shells 


Ccordenser: 


Velocity 


Eressure 


Jube length 


Shell dia. 


lube dia. 


# of skells 


Mass flow rate 


System ccst 


Range 


J. 32066 Soe 
(1.013-1.029) 
aferox Q 


39. 16-40 .00 
(11.91-12.19) 
24.12-25.00 
(7. 352-7.620) 
0.756-0.761 
(1.919-1. 934) 
7.141-7. 409 


2.793-2.823 
(0. 851-0. 860) 
91.53-91.60 
(631. 1-631.6) 
39-35-40 .00 
(11.99-12.19) 
44 665-45 00 
(13 .6 2-13.72) 
0.7€32-0.7686 
(1. 938-1. 952) 
3.297-3.401 
gz 1e= 22001 
(10. 03-10.07) 
7.57&-7.673 


Table 6 
EXAMEFLE PHOBLEYA SOLUTION SUMMARY 


Icx 


ean 


3.335 
(1.016) 
119.2 
(815.0) 
39.71 
(12.10) 
24.91 
(7.593) 
0.758 
(1.926) 
7.228 


2.808 
(0.856) 
91.56 
(631. 3) 
39.78 
(12.12) 
44.89 
(13.68) 
0.7660 
(1.946) 
a4 
227 
(10.06) 
7.607 
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Vax 
x 


Jnits 


ft/sec 
(m/sec) 


(Ca) 


ft/sec 

(n/SEc) 

lbyin? 

(Kp) 

Ew 

(m) 

fe 

(0) 

in 

(cm) 
10¢lbmyhr 

(10%kg/kr) 
106 $ 





The arithmetic mean of the seven soluticrs is calculated. 
The difference between the maximum and the minimum valves is 


shown as a fercent of the mean. The percent variaticn gives 


sche idea cf how well the solution was determined. Note 
that the tuke length and shell diameter constraints are 
binding for ktoth the boiler and the ccndenser. The 


sukprceblem still terminates only when the penalty furction 
does not decrease from one iteration to the next. The 
Magritude cf Wy P (X,T) iS not becoming small enough to 
satisfy tke ccnvergence critericn. The value never beccmes 
less thar C¢C.02 and is usually greater than 0.5 when a 
sukgreblea is stopred by the progran. 


Useful results are obtainable, however. The encineer 
can learn much about his problem even from a single _ foint. 
Fer exagfple, ‘frcem the cost breakdown in Figure &€, the 
portion cf tke costs attributable to each corponent can be 
seen. The heat exchangers, of course, dominate with 78% of 
the cost. Tke circulation pumps are about 21% of the cost 
with the feed pumps being 1.2% of the total cost. Since 
the efficiency assumed for the circulation fumps is frebably 
toc high, these pumps would be even mcre ccstly in terms of 
energy and mcney. This could indicate the priorities that 
should ke set on the design efrorts for improvements cf the 
Varicus ccmpenents with the highest priority beinc the 
highest percentage cost component. 


The systents fewer distribution 1s shown in fable 7. 
The sigrificant features are the large energy akscrrtion 
rate for the saall useful output fower and the parasitic 
rusping fower. The lower the efficiency the greater must be 
the power aksorbed for the same output. This is to be 
expected fer this kind of system. The parasitic priamfing 
rewer (12% of the turkine output) is aktout twice what it is 


for a fcessile fuel plant. 
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CONDENSER 


BOILER 
$3,333,000 $2,528,000 
43.9% 33.3% 


CONDENSER 
CIRCULATION 
PUMPS 
$978,000 





PRED BOILER 


PUMPS CIRCULATION 12.9% 
$94,100 PUMPS 
ie. 25 $654,000 

8.6% 


FIGURE 8. COST BREAKDOWN 
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Table 7 
SOLUTICN EOWER BUCGET 


BIUZbr MWatts 
Beat Aksorption 3.56x 109 1042.8 
Turbine Shaft Power 9.69x107 28 24. 
Feed Fumps 4.49x106 Tae 
Elxr Circulaticn Pumps 2.58x 106 0.8 
Ccnd Circulation Pumps 4.51x106 i 
Feat Bejected | 3.47x1C9 1901527 
Net Useful Power Bs 54010? 2 gest) 


68 





Cnce the mcdel has been optimized with the number of 
shells as real variables, the problem can Le re-solved twice 
by fixing tke shells at the next higher and the next lower 
integer numker. For the higher number of shells, the cost 
increased slightly as expected, from 7.607 million dcllars 
average to 7.983 million dcllars. The problem with the 
lower integer number of shells could not be sclved dve to 
the constraints cn the tube length and the diameter cf the 
shells keing binding before the shell numbers were 
integerized. The program is unable to find a new solution 
that satisfies the equality constraints. 


Table 8 shows some miscellaneous data from the mcdel 
that deserves comment. The overall heat transfer 
ccefficient cf the boiler is approximately three times that 
for the ccndenser which is the major reason for the lcwer 
ccest of tke boiler relative to the condenser. However, 
rememrer that the bundle factor, egual to 3, is a part of 
the heat transfer coefficient. More data must ke gotten 
akcut this factor in crder to improve the reliance on any 
sclution. There is a significant difference between the 
Prcpane heat transfer coefficients that were assumed ard the 
theoretical values calculated from the correlations. Other 
values shculd be tried until the fixed value is withir the 
20% tolerance on the theoretical value. The flcw of 
seawater thrcugh the tubes in both the kciler ané the 
ccrdenser is turbulent. Apparently, the resultinc cost 
Savings due to the improved heat transfer rate through a 
smaller area is more than the increased costs of the higher 
Fumping power needed. It seems likely that the solution 
oktained is globally optimal; several different starting 
Fceints were tried. The problem has been started in Eoth the 
lawinar and the turbulent flow regions with no change in the 


results. 
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Tabie 8 


ADDITIONAL CATA FROM SAMELE 
PROBLEM SOLUTION 


Working Fluid Heat Transfer Ccefficients 


Assuned 145. 
(0 2182 2) 
Theoretical 227. 
CV229) 
Surface Area per Shell 2239005 
(2.22x10*%) 
Seawater Cutliet Temp eis. 3 
(22.4) 
Shell Saturaticn Temp 66.7 
CVS) 
Seawater Mass Flow Rate 4.39x108 
(1.99x108) 
Mean Temp Difference 6.67 
(3.71) 
Overall Heat Transfer 310.6 
Coefficiert (1.76) 
Heat Flux 2070.0 
(6.520) 
Seawater Keyrclds ro. 19,800. 
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145. BTU/hr-ft2-9P 
(0.822) (kW /m2-°C) 
241. 
(1.37) 
1.27x106 £t2 
(1.18x105) (m@) 
43.8 oP 
(6.55) (°C) 
49.9 Sis 
(9.94) (°C) 
9.11x109 lpom/hr 
(4.13x109) (kg/hr) 
7.85 OF 
(4.36) (°C) 
105.0 ETU/hr-ft2-9F 
(0.595) (kW¥/m--9C) 
824.9 BQU7 feo -ic 
(2.602) (kW/m2) 
97960. - 





IV. RECOMMENDATICNS 


SE ee a ae eee ee ce 


The present model has a great deal mcre informzaticn to 
offer, but scme improvements in the model shculd be made te 
improve the ability of a nenlinear program to solve the 
prccblem. The first recommendaticn would be to study the 
equality constraints and place inequality ccnstraints cn any 
variable cr grouping of variables if they Ecth can tkecome 
negative and also must have their logarithms comguted. 
Examples, already included in the program, are the numter of 
tukes in tke shells. The heat transfer coefficients fcr the 
working fluid side in the boiler and the ccndenser cculd be 
included as eéquality constraints. This can be dcne by 
adding the heat transfer coefficients to the vector xX and 
réeguiring that the variables be equal to the coefficients as 
ccmputed Ly the appropriate correlation. The same cculd be 
dcne with FEalen and Small's property-dependent calculation 
of the sefaraticn area by equation (12) and (13). Frcm an 
cverall standpoint, the results seem to indicate either the 
Ereblem is tcc difficult to be solved efficiently ir its 
present form or that another ncnlinear prcegram, other than 
SUMT, might te better suited to the problen. Numerical 
differentiaticn of the constraints in this froblem de not 
seem to help the program to find the solution since the 
magnitude cf WP(xX,r) is net near zero when a sukfroklen 


SWerTs. 


Once the optimization scheme is functioning froferly, 
the model can te used to study the effects cf fouling and 
cerrosicr ané tne costs of their prevention. Many profosals 
have teen made that must be analyzed in the ccntext of the 
total systen. Scme such aS Amertap, where plastic fcam 
kalls are pumped through the tubes to clean the heat 
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transfer surfaces, adda cost penalty fron the capital and 
cperating ccst of the Amertap system and from the increased 
Fumping fCwer required. Coatings add the capital cests of 
the coatirg and the added capital cost of increased heat 
transfer area required hecause of the increased tkermal 
Tesistance cf the ccating. 


A sensitivity analysis can te performed on the model to 
test the assumpticns and the various ccnstants. The 
sensitivity te changes in these factors can indicate how 
much one shculd be willing to pay for the improvement 
sought. High sensitivity indicates areas where the greatest 
effort for inprovement should Fe desired or the priorities 
Shculd Ee set cn research into the unknewn areas cf the 
problem. The trade-off between cost and performance can be 
studied Eby changing the output power and looking fcr the 
lowest ccst fer kilowatt-hour, which is the measure by which 
corventicnal fceower plants aré compared. One could icok into 
the cost and the Fenefits of increasing the depth from thich 
the cold water is pumped versus the improved temperature 
differential fpossiktle from usirg the colder water. 


For the capital cost of the system to be complete, the 
turkine design must be included and ccst e¢stimatés made. 
The turktine will be unique due to the very low pressure 
rotential drcp it must use. Some work has been dore by 
cther grcups working on OTEC, Ref. 25 and 26. Reference 27, 
by Rcbert L. Bartlett, can provide costing and design 
infcrmaticn although he deals only with steam turbines. The 
cremaining fart of the system that is desireable FOr 
inclusicr ir the analysis presented heréin is the cold 
seawater intake pipe. It will be on the order of 1,0C0O it 
(304.&m) lcrg and 100 ft (30.4 m) in diameter. 
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